Signature of a universal statistical description for drift-wave plasma turbulence 
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This Letter provides a theoretical interpretation of numerically generated probability density func- 
tions (PDFs) of intermittent plasma transport events. Specifically, nonlinear gyrokinetic simulations 
of ion-temperature-gradient turbulence produce time series of heat flux which exhibit manifestly non- 
Gaussian PDFs with enhanced tails. It is demonstrated that, after the removal of autocorrelations, 
the numerical PDFs can be matched with predictions from a fluid theoretical setup, based on the 
instanton method. This result points to a universality in the modeling of intermittent stochastic 
process, offering predictive capability. 



Stochastic physical processes are most often observed 
to be unimodal with exponential tails [l[ , a feature which 
is also attributed to fluctuations in magnetically con- 
fined plasmas Q- Q. These fluctuations are intermittent 
events manifesting a patchy spatial and bursty tempo- 
ral structure, pertaining to radially propagating coherent 
structures like blobs or avaloids [8|, and have been sug- 
gested to carry a significant fraction of the total trans- 
port Q . Therefore, a comprehensive predictive theory is 
called for, in order to understand and, subsequently, im- 
prove properties related to intermittency. A major goal 
would be, for instance, the control of edge heat flux loads, 
which depend on the instant amplitude of fluctuations, 
as opposed to the mean load, which can be calculated by 
quasilinear theory. 

In terms of mathematical description, the likelihood of 
intermittent events related to plasma turbulence is ex- 
pressed by probability density functions (PDFs), which 
usually deviate significantly from the Gaussian distribu- 
tion. Along these lines, there have been attempts to char- 
acterize the statistical properties of the PDFs, based on 
phenomenological premises [13, El or numerical investi- 
gations [H, [l3[ alone. Here, however, we carry out a di- 
rect comparison between first-principles analytical mod- 
eling and numerical simulations. As will be evident in 
the sequel, although the two approaches express the same 
physics, they nevertheless greatly differ in their theoret- 
ical backgrounds. A key finding of this work is that the 
intermittent process in the context of drift-wave turbu- 
lence appears to be independent of the specific model- 
ing framework, opening the way to the prediction of its 
salient features. 

The main part of this Letter consists in providing a 
theoretical interpretation of the PDFs of radial heat flux 
derived by nonlinear, local, gyrokinetic (gk) simulations 
of drift wave turbulence in tokamaks. Our paradigm 
in this work is Ion- Temperature-Gradient (ITG) turbu- 
lence with adiabatic electrons [l4| ■ The simulations have 
been carried out with the GENE code [lj| in a simple 
large aspect ratio, circular tokamak geometry. In par- 
ticular, we calculate the turbulent ion radial heat flux, 



f v r J (h m i v '\\ + l 1 ^ J Fij , where v r is the radial 
E x B velocity, v» the parallel velocity, \i the magnetic 



moment, rrii the ion mass, Fi the perturbed ion gyrocen- 
ter distribution function and B the modulus of the mag- 
netic field. The brackets denote spatial averaging over 
the entire simulation domain. In order to perform a re- 
liable comparison, we produce a series of different cases, 
by varying the magnetic shear s — (q is the safety 
factor, i.e., the number of toroidal turns of the magnetic 
field for each polodal turn), which appears both in the 
gk simulations and the theoretical model. There is noth- 
ing special about the selection of this parameter, and a 
similar scan could be performed over, for example, the 
ion temperature or density gradient. Here, we have set 
the (normalized) ion temperature gradient R/L^ = 9, 
the density gradient R/L n — 2, and the ion-to-electron 
temperature ratio r = 1. 

For each case, the time evolution of the ion heat flux 
is considered as a time series, to which we apply stan- 
dard Box-Jenkins modeling 16]. This mathematical pro- 
cedure effectively removes deterministic autocorrelations 
from the system, allowing for the statistical interpreta- 
tion of the residual part, which a posteriori turns out to 
be relevant for comparison with the analytical theory. In 
our setup, it turns out that an ARIMA(3,1,0) model ac- 
curately describes the stochastic procedure, in that, one 
can express the (differenced) heat flux time trace in the 
form 

Qt+i = ai Q t + a 2 Qt-i + a 3 Q t _ 2 + Qres(t) (1) 

where the fitted coefficients 01,02,03 decribe the deter- 
ministic component and w(t) is the residual part (noise). 
The reliability of the produced autoregressive model was 
routinely verified via portmanteau testing and ovcrfit- 
ting. 

It is systematically observed that the residual PDFs 
are manifestly non-Gaussian with elevated tails, as shown 
for instance in Fig. [I] for the case s = 1.0, where the sam- 
ple residuals from a GENE simulation are tested against 
the Gaussian distribution, via a normal quantile-quantile 
plot. Furthermore, two significant statistical quantities, 
namely the variance a 2 — (Q 2 ) — (Q) 2 and the kurto- 
sis K — (Q 4 ) /c 4 (here, the brackets denote averaging 
over the statistical sample), during the magnetic shear 
scan are summarized in Table 1. Evidently, the decrease 
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of the variance and the rapid increase in kurtosis for in- 
creasing magnetic shear renders a Gaussian description 
improper. We note in passing, that the values of s were 
selected such that the resolution of the simulation domain 
remains constant, thus avoiding numerical artifacts. 

In the sequel, it will be shown that the aforementioned 
PDF tails can be satisfactorily predicted, starting from 
a fluid model consisting of a continuity and an energy 
equation (l7j . and using a nonperturbative statistical 
technique, called the instanton method, which has been 
adopted from Quantum Field Theory and then modi- 
fied to classical statistical physics for Burgers turbulence 
and in the passive scalar model of Kraichnan [TJ[. In 
this context, the key element is to identify the bursty 
or intermittent event with the appearance of a coherent 
structure (e.g., streamer). 
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FIG. 1: Normal quantile-quantile plot showing strong devia- 
tion of the residual PDF from Gaussian statistics, in view of 
the enhanced tails representing about half of the sample. 



TABLE I: Simulation Results 



s Variance Kurtosis 

0.25 5.962e-4 3.529 

0.40 5.253e-4 

0.50 6.049e-4 

0.60 5.874e-4 

0.75 4.403e-4 

1.00 8.253e-5 



3.517 
3.527 
3.715 
4.324 
14.413 



In the following paragraph, we briefly outline the im- 
plementation of the instanton method. For more details, 
the reader is referred to the existing literature [2(|- [23| . 
The PDF tail is first formally expressed in terms of a path 
integral by utilizing the Gaussian statistics of the forcing 
in the continuity equation in a similar spirit as in Q . An 
optimum path will then be associated with the creation of 
a modon (among all possible paths or functional values) 
and the action (S\, below) is evaluated using the saddle- 
point method on the effective action in the limit A — > oo. 



The instanton is localized in time, existing during the for- 
mation of the modon. The saddle-point solution of the 
dynamical variable cj>(x, t) of the form <fi(x, t) — F(t)ip(x) 
is called an instanton if F(t) = at t = — oo and F(t) ^ 
at t = 0. Note that, the function ip(x) here represents 
the spatial form of the coherent structure. Thus, the in- 
termittent character of the transport consisting of bursty 
events can be described by the creation of modons. The 
probability density function of the heat flux Q can be 
defined as 

P(Q) = (J(« r nT i (f= x )) -Q))= [ d\e lXQ h, (2) 



where 



h = (exp(-iXv r nTi(x = x ))). 



(3) 



Here, v r is the radial drift velocity and Ti the ion temper- 
ature. The integrand can then be rewritten in the form 
of a path-integral as 



V4>V<t>e 



(4) 



Although <f> appears to be simply a convenient mathe- 
matical tool, it does have a useful physical meaning that 
should be noted; it arises from the uncertainty in the 
value of due to the stochastic forcing. That is, the 
dynamical system with a stochastic forcing should be ex- 
tended to a larger space involving this conjugate variable, 
whereby and constitute an uncertainty relation. Fur- 
thermore, acts as a mediator between the observables 
(heat flux) and instantons (physical variables) through 
stochastic forcing. In Eq. [4j the integral in A is com- 
puted using the saddle-point method where it is shown 
that the limit A — > oo corresponds to Q — » oo, repre- 
senting the tail part of the distribution. Based on the 
assumption that the total PDF can be characterized by 
an exponential form, the expression 
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is found [2J], where the heat flux Q plays the role 
of the stochastic variable, with P(Q) determining its 
statistical properties. Several auxiliary definitions are 
also utilized; the normalization constant N; the gra- 
dient scale lengths Lf = — (dlnf /dr) 1 ; the normal- 
ized modon speed U = RU/L n and temperature ratio 
t = Ti/T e ; g l = Lj d /uj± = 2^l(cos(£) + s£sin(£)) where 
u!d is the curvature drift frequency and is the dia- 
magnetic drift frequency; k\ = kz{l + s 2 £ 2 ) is the per- 
pendicular wave number; the brackets denote averaging 
along the field line, e.g. for an arbitrary scalar func- 
tion /, (/) = \\ <£/</>/ W d£0 2 where the eigenfunc- 
tions 0(£) are extracted from the GENE simulations and 
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A 4 = (Q) i s the mean value of the heat flux. The coeffi- 
cient 60 is a free parameter and represents the strength 
of the forcing in the continuity equation. Note that the 
proposed PDF is close enough to a Gaussian distribu- 
tion to match the bulk of the PDF while retaining the 
fat tails. Furthermore, the exponential form of the PDF 
will be the same for momentum flux with a modified co- 
efficient b, which is in agreement with findings in recent 
experiments [24| . 

In Fig.H(i = 0.25), the PDF of heat flux Q stemming 
from the raw timetrace of the simulation (blue line) is 
shown. Although a simple visual inspection precludes 
the Gaussian distribution, the analytically predicted in- 
stanton method from Eq. [5] (red line) does a rather poor 
job modeling the simulation result. This is, however, 
contrasted to the PDF using the weak non-linear model 
(dashed-dotted black line) where reasonably good agree- 
ment is found similar to what have been reported ear- 
lier (y, |9j. The weak non-linear model (a Laplace dis- 
tribution) is also predicted by the instanton model [25[ 
when the non-linearities are neglected and has the form 
P(Q) = e~' < 5 _A1 '/ b /(2fo) where b is determined by the vari- 
ance. The coefficient bo in Eq. [5] is determined from the 
simulations at the point s = 0.5. 



between the gyrokinetic and fluid descriptions of drift 
wave turbulence strongly suggests a universality of the 
statistics of the intermittent process. 
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FIG. 3: Comparison of the numerically estimated PDF 
(blue line) against the analytically predicted instanton model 
(dashed-dotted black line) and weak non-linear model (dashed 
green line) compared to the best fit Gaussian (dashed red line) 
for s = 0.25. 
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FIG. 2: The numerically estimated PDF of heat flux (dashed 
blue line) in comparison with the analytically predicted in- 
stanton method (red line) and the weak non-linear model 
(dashed-dotted black line) at magnetic shear I = 0.25. 

Figs. [3] and S] single out the comparison of the numeri- 
cally estimated residual PDFs against the instanton pre- 
diction in Fig. [T] Here, the PDF of Qres from the sim- 
ulation (blue line) and the analytically predicted PDF 
from the instanton method (dashed-dotted black line) 
demonstrate a significantly better agreement both at the 
tails and the center of the PDF, as compared to the best 
fit (using the first two moments) Gaussian distribution 
(dashed red line) at magnetic shears s = 0.25 and s = 1.0, 
respectively. Note, we show the weak non-linear model 
with a value of the variance of 8.00e-4 for better agree- 
ment and that the distribution does not capture the form 
of the stochastic PDF. The surprisingly good agreement 
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FIG. 4: Comparison of the numerically estimated PDF (blue 
line) against the analytically prediction (dashed dotted black 
line) and the best fit Gaussian (dashed red line) for s = 1.0. 

Since the quality of the agreement in the previous fig- 
ures relies heavily on the exponential coefficient b in Eq. 
5, this is displayed for the whole magnetic shear scan in 
Fig. [5J verifying that good agreement is expected for any 
value of s. Here, we note that the abrupt bending of 
the curve for larger values of shear is solely attributed to 
FLR effects. 

In conclusion, we have presented a first quantitative 
paradigm of universality in intermittent stochastic pro- 
cesses related to drift wave turbulence. Numerical PDFs 
of heat flux were generated with the gyrokinetic code 
GENE in the framework of toroidal ITG turbulence, and 
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FIG. 5: The coefficient b as a function of magnetics shear 
s. The simulation results (red squares) are compared to the 
analytical fluid model prediction (black stars) and the same 
with suppressed FLR effects (blue diamonds). 



firm the exponential form of the PDFs, therefore adding 
the important element of predictive strength to the ex- 
isting phenomenological approaches. 

In future publications we will address the emergent 
universal scalings of the PDFs of potential, density and 
temperature where the theory predicts different tails [25[ . 
A study of potential or density fluctuations opens up op- 
portunities for comparison with experimentally measured 
PDFs. Still within the scope of ITG turbulence, we will 
include effects from kinetic electrons, in order to test the 
robustness of the exponential scaling. A quite interesting 
complementary work would involve incorporating results 
from nonlinear fluid simulations, based on Braginskii-like 
equations, and relate them to the already presented find- 
ings. Finally, our present setup provides a testbed for the 
investigation of the impact of zonal flows on the statis- 
tics, as reported elsewhere [Hj]. Work addressing these 
avenues has been initiated. 



subsequently processed with Box- Jenkins modeling, in 
order to remove deterministic autocorrelations, thus re- 
taining their stochastic parts only. These PDFs have 
been shown to agree very well with analytical predic- 
tions based on a fluid model, on applying the instanton 
method. Specifically, we were able to quantitatively con- 
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